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CUSP EXCURSIONS FOR THE EARTHQUAKE FLOW ON THE 
ONCE-PUNCTURED TORUS 


SER-WEI FU 


Abstract. In this paper we study the typical speed of a generic earthquake 
trajectory leaving compact sets in the moduli space of the once-punctured 
torus. Mirzakhani showed that the earthquake flow is measurably equivalent 
to the horocyclic flow, which has been studied extensively. Our main result 
shows that the earthquake flow and the horocyclic flow behave very differently 
in cusp excursions. In particular, we prove a relation between the systole 
function and continued fractions and discuss the cusp excursions of earthquake 
trajectories. 


1. Introduction 


The motivation of the paper comes from the result of Mirzakhani in [8] showing 
that the earthquake flow is ergodic with respect to the natural finite measure v(S) 
coming from the Weil-Petersson measure on the bundle P 1 A i(S) of unit length 
geodesic measured laminations. Throughout the paper we will let S = S\,± be the 
once-punctured torus, A i(S) be the moduli space, and T(S) be the Teichmiiller 
space. We use E t (X. A) to denote the earthquake flow and i sys (X, A) to be the 
length of the systole of X (See Section [2] for definitions). 

The study of cusp excursions for the earthquake flow is aimed at understanding 
the constant 


Cf(X , A) = limsup 

t—>oo 


-log e sya (E t (X, A)) 

m 


Ergodicity implies that Cf(X, A) is a constant for almost every ([A], A) with respect 
to v(S). The goal is to find the threshold function / such that Cf is a nonzero 
number that depends only on the surface. 

One of the classical approaches to cusp excursion is the logarithmic laws initiated 
by Sullivan m- The Borel-Cantelli lemma and a volume estimate of the thin part 
of the moduli space can be used to obtain a discrete upper bound for C\ og n\ . More 
specifically, if the measure of the e-thin part is given by g(e ), then we obtain 


lim sup 

n—> oo 


-\og£ sys (E n (X,X)) 
log n 


< lim 

£—>-0 


lQ g(g(g)) 

loge 


where n G N, 


for almost every ([A'], A) G V 1 M.(S g ) with respect to v(S), see Sectionjd] 

The logarithm law of the horocyclic flow proved by Athreya in JT] states that 
the discrete upper bound above is actually an equality even in the continuous case. 
However, the earthquake flow fails to satisfy the same property. Our main result 
shows that is infinity almost everywhere for the once-punctured torus. 


l 
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Theorem 1 . Let S be the once-punctured torus. For almost every ([A], A) £ 
V 1 A4(S) with respect to v(S), we have 


lim sup 

t—>oo 


-log e sys (E t (X, A)) 

Vt 


= 00 . 


The theorem echoes with the result shown by Minsky and Weiss [7, Proposi¬ 
tion 8.1] stating that in the case of the once-punctured torus, the trajectories of the 
horocyclic flow and the earthquake flow can be infinitely far apart. As a corollary 
of the proof of the main theorem, we identify a family of earthquake trajectories 
associated to the golden ratio that goes arbitrarily deep into the thin part with the 
slowest possible speed. 

We study the systole function along earthquake trajectories by explicitly parametriz¬ 
ing the space of hyperbolic metrics, the simple closed geodesics, and the space of 
measured laminations on the once-punctured torus. The systole functions on earth¬ 
quake trajectories are closely related to the continued fraction expansion of slopes 
that depends on ([A'], A), see Section [3] We show the result of the classical argu¬ 
ment for cusp excursions in Section [4] for completeness and prove the main theorem 
in Section [5] The proof is a straight forward computation of the limit supremum 
for any given (A, A) and the measure of the set that satisfies our statement. 
Acknowledgement. The author thanks Jayadev Athreya for suggesting the project 
as well as many useful conversations. The author is also grateful for discussions 
with Christopher Leininger, Matthew Stover, and Joseph Vandeliey. 


2. Background 

We use a similar choice of notation as in [8]. Let T{S) be the Teichmiiller space, 
which is the space of hyperbolic metrics on S. The mapping class group Mod(S'), 
the group of automorphisms of S up to isotopy, acts on T{S) and the quotient 
is the moduli space A4(S). The Weil-Petersson symplectic form induces a finite 
measure on A i(S). 

Let A4C(S) denote the space of measured laminations on S (see [2] ). The space 
A ijO(S) carries a mapping class group invariant volume form fiTh introduced by 
Thurston. Given X £ T(S) and A £ AiC(S), we denote the length of the lamination 
A in the hyperbolic metric A' by £\(X). Finally we let VT(S) = T{S) x MC{S ) be 
the bundle of geodesic measured lamination over T{S), V 1 T(S) be the unit sub¬ 
bundle with respect to the hyperbolic length of A, and V 1 M(S) be the quotient 
V 1 T(S)/ Mod(S). The earthquake-invariant measure on V 1 M(S) considered by 
Mirzakhani in [8] is a finite measure denoted by v(S). 

The earthquake map is a generalization of twisting along simple closed geodesics 
with respect to a hyperbolic metric on S , see Emu. We denote the earthquake 
flow at time t by E t (X, A) £ VT(S), where E t ( A, A) = (Y, A) with Y being the 
hyperbolic metric obtained by twisting A along A for distance t. 

The earthquake flow on V 1 Ad(S) is ergodic with respect to v{S). This implies 
that almost every trajectory of the earthquake flow in Ad (S) is dense. Alternatively, 
we consider the systole length function £ sys that maps (A, A) £ VT(S) to the length 
of the shortest closed geodesic with respect to A. The ergodic property implies that 

lim inf £ sys (E t ( A, A)) = 0 for almost every (A, A) £ 'PT(S). 
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The Teichmiiller space of the once-punctured torus T(Si,i) is isometric to H 2 
and is parametrized via the Fenchel-Nielsen coordinate. We fix a pair of generators 
for the fundamental group of S ip that intersect once. Then every oriented simple 
closed geodesic corresponds to a pair of integers given by the algebraic intersection 
number with the generators. We denote an oriented simple closed geodesic by 
7(0,6), where gcd(a, 6) = 1. 

The space of measured laminations is parametrized by the slope (the 

ratio of a and 6) and the hyperbolic length with respect to X. We say that a mea¬ 
sured lamination is rational if it is a weighted simple closed geodesic and irrational 
otherwise. 

Let us setup the Fenchel-Nielsen coordinate. For X £ T{S 17) and 7 = 7(0,6) 
a simple closed geodesic, the Fenchel-Nielsen coordinate of X with respect to 7 is 
(£ 7 (AT), t 7 (X)), where r 7 (X) is the twist parameter. We choose the twist parameter 
in the following manner. Let 7(c, d ) be an oriented simple closed geodesic such that 
ad — 6c = 1 . The twist parameter r 7 (X) = 0 if and only if X is the hyperbolic 
metric that minimizes the length of 7(c, d) for the fixed length £^(X ). The twist 
parameter is normalized such that integer valued twist parameters correspond to 
metrics obtained by Delm twists about 7. The twist parameter r 7 increases linearly 
along the earthquake trajectory of E t (X, 7). 

3. The punctured torus 

In this section we prove that the systole length function on a earthquake trajec¬ 
tory can be approximated via continued fraction expansions of irrational slopes. 

Proposition 2. Fix a simple closed geodesic 7. If X £ 7~(Si,i) has Fenchel-Nielsen 
coordinate ( L,p/q ) with respect to 7, where gcd (p,q) = 1 and p/q £ (0,1), then 

C x {L)e- L l ( 2 ®> < £ sys (X) < C 2 {L)e~ L /^\ 

where C\(L), C 2 (L) both limits to 4 when p, q are fixed and L goes to 00 . 

Proof. We consider a once-punctured torus with a hyperbolic metric A' £ T(5'i,i). 
Figure [I] below shows a fundamental domain in the upper half plane that is chosen 
canonically with a given Fenchel-Nielsen coordinate (L,p/q) with respect to 7. 



The edges are identified by colors via hyperbolic isometries. The length ratio 
between the blue and the orange edges is equal to p : q — p. We use R{L) to denote 
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the length of the shortest geodesic arc with endpoints on 7. It is well known that 
(R o R)(L) = L and 

/ e £/2 + i \ 

i?(L) = 2 1 °g(^- Z 7^± T j. 

For fixed p/q £ ( 0 , 1 ), the systole function is realized by a simple closed curve 
that intersects 7 more than once when L = £^(X) is large enough. We use the 
analogy from the fiat torus case and look at the length of simple closed curves a 
that intersect 7 exactly q times as seen in Figure [ 2 ] 


oi\ 



_ 1 _ 1 _ ; _ ; _ 


Figure 2 . The flat torus analogy. 


We obtain the following inequality involving a, the length of a, with respect to 
L which came from switching the role of 7 and a in the pictures. 

qR(a) < L < qR(a) + 


Since 7 consists of q geodesic segments from a to a, the left-hand-side is obtained 
from R(a). The riglit-liand-side is obtained by considering a piecewise geodesic 
homotopic to 7. The piecewise geodesic goes along the geodesic segment that 
realizes the injectivity radius of a q times and travel along a otherwise. 

Straight forward computation shows that 


=-“/ 4 


1 + 


< e- L '^ < 




1 + 




A closed curve homotopic to a is obtained by concatenating q radial segments 
between the two circles in Figure Q] The hyperbolic length of the radial segments 
goes to zero as L —> 00 since the Euclidean length of the segments goes to zero 
while being bounded distance away from the x-axis. Hence a — > 0 as L —► 00, 
which implies that e -a / 4 is approximately 1 and e a / 2 — 1 is approximately a/ 2 . In 
other words, there exists C\(L ), that both converges to 4 such that 

C'i(L)e- i/(2<?) < a < C 2 (L)e~ L ^ 2q) . 


To complete the proof, we need to show that the length of a is shorter than any 
other simple closed geodesic for large L. In fact, the geodesic homotopic to a is the 
only simple closed geodesic whose length goes to zero as L —> 00. Therefore the 
estimate holds for the systole length. □ 


The estimate given by Proposition [ 2 ] allows us to construct explicit sequence of 
sample points along the earthquake trajectory. The sequence of sample points will 
be used in the proof of the main theorem for a lower bound of C^(X, A). 

Proposition 3. For any X £ T( 5 i j i) and any irrational A £ AiC(S\ } \), we have 

lim inf i sys E t ( X , A) = 0 . 
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Proof. It suffices to construct a sequence {t n } such that £ sys Et n (X, A) goes to zero 
as n —> oo. Without loss of generality, we assume that £\{X) = 1. There exists a 
sequence of simple closed geodesics {7 n } such that (T Tn (A)) _1 7 n —> A in At£(Si,i). 
Fix o > 0 and T > 0. There exists N such that n> N implies the following. 



Let a be a closed geodesic that realizes £ sys {X, A) and 7 = 7 m for some m> N. 
We will assume that m is large enough from here on, hence the length of 7 will 
also be large enough. We observe that there exists po,qo such that a is the systole 
of Xo = ( £-y(X),po/qo ) in Fenchel-Nielsen coordinate with respect to 7. Let X\ 
be (£j(X),pi/qi) with pi/qi be the least rational number larger than po/qo with 
1 < qi < qo. We iteratively define Xj from Xj_\ similarly and choose tj to 
correspond with Xj on the earthquake flow. Let k be such that tfe < T < tk+i- 
The choice of Xj’s along with the estimates in Proposition[2]imply that the finite 
sequence 


{IsysE^X^iX))- 1 ^ 


is decreasing. Therefore 


inf £ sys E t (X,X) < (4 + e)e-^ (A ' )/(2 « fe) 


0 < t<T 


The value of qu uniquely depends on po, go, and k. Since the length of a is fixed, 
we conclude that the ratio £ 1 {X)/qo as a function of m converges to a constant. 
Hence it suffices to consider the ratio go /qk and show that it goes to infinity as m 
goes to infinity. 

We use some properties of continued fractions from [4] in the rest of the proof. 
By construction, for qj ^ 1, Pj/qj is a convergent of po/qo- To be more specific, if 
Po/qo is [ao; ai, 02, ■ ■ ■, a n \ in the usual convention of continued fraction expansion, 
then pj/qj is [do; 01,02,..., a^Q j], where h(j) is equal to n — 2 j or n — 2j + 1 since 
only the odd convergents are greater than po/qo- 

The ratio qo/qk depends on the continued fraction expansion of po/qo ■ In partic¬ 
ular, the ratio grows with k unless qk = 1. We observe that qo/qk goes to infinity 
as k goes to infinity whenever po/qo converges to an irrational number as m in¬ 
creases, which is a consequence of A being irrational. Therefore for any fixed e, as 
T increases, m and k go to infinity and that completes the proof. □ 


4. Discrete upper bound 


We include a discussion on the classical upper bound obtained by the Borel- 
Cantelli lemma for completion. 

Proposition 4. Let be the finite measure on V 1 A 4 (S ip). Then 



Proof. Recall the McShane identity and the computation from [9] showing how the 
volume of At 1,1 is computed using Theorem [5] in [6]. 
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Theorem 5 (McShane identity). Let X be a point in the Teichmuller space of the 
once-punctured torus. Then we have 



7 


where £ 7 (X) is the length of the geodesic 7 with respect to X and the sum is over 
all simple closed geodesics 7 on X. 


r 00 f°° £ 7 T 2 

Vol(X ltl ) = 2 / if(l) dt= 2 —- dl=—. 

Jo Jo 1 + e D 

In [8] Mirzakhani showed that the earthquake-invariant measure v restricted to 
the space of unit-length measured laminations is equal to the Thurston volume of 
the unit ball in the space of measured laminations. Therefore we can compute the 
volume of the £-thin part of V 1 with the help of the following claim. 

Claim. There exist constants C\,Ci > 0 such that for any X G T(Si,i) with 

^sys(X ) = £, 

-r— < HTh{\ e MC(S lA ) I h(X) < 1} < . 

£ log £ £ log £ 

Proof of Claim. Let 7 be a systole of X. The Fenchel-Nielsen coordinate of X with 
respect to 7 is (£, r) for some r. Let a be a simple closed curve that intersects 7 
once. Then A £ MC{S) has coordinate (i(A, 7), i( A, a)) where t(-, •) is the algebraic 
intersection number with some fixed orientation. 

The Thurston volume can be computed by a double integral of the characteristic 
function of the unit ball in Since the length of A is linear with respect 

to its coordinate, we can rewrite the integral in polar coordinates. 

ft pr(6) 

I»tk {A G M£(S) I h(X) < 1 } = / / rdrdd , 

Jo Jo 

where r( 9 ) is the reciprocal of the length of the measured lamination with coordinate 
(cos 9 , sin 6 ). 

The constant bounds controls the following estimates. We can choose a such that 
J a (X) is bounded between R(e) and R(e) +£, where R(e) is as defined in the proof 
of Proposition [ 2 ] The length of the measured lamination (cos 0 , sin 6 ) is bounded 
between max{I?(£)| cos 0 |,£| sin 0 |} and (/ 2 (e)| cos 0 | + e\ sind|). The integral 

r 1 

Jo (i?(£)| COS0| + £| sind|) 2 ^ 

is estimated by the standard method of breaking it into pieces at the points when 
i?(£)| cos 0 | = e\ sind|. Finally we estimate R(e) by log(l/£) when £ is small. 

That finishes the proof of the claim. Finally, we consider the integral below and 
standard computation will complete the proof. 

/ ' TT ~1 d£ = O ( ] as £ -)• 0. 

Jo 1 + e e £ log t \\ogs) 

□ 

The following proposition is the result of a standard argument using the Borel- 
Cantelli lemma along with the volume of the thin part as stated in Proposition | 4 ] 
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Proposition 6. Let S = S i,i and t{n) be any fixed sequence of real numbers. For 
almost every (X, A) £ T(S) x XiX(S), we have 


lim sup 

n—>oo 


-log £sys(E t (n) {X, A)) 
logn 


< 1 , 


where £ sys is the systole function and Et is the earthquake function for time t. 


Proof. Recall the easier part of the Borel-Cantelli lemma. 


Theorem 7 (Borel-Cantelli lemma). Let (X, X, p) be a probability space, where 
X is the sample space, T is a set of events, and p is a probability measure. Let 
X n : X —> {0, 1} be a sequence of random variables with p{x £ X : X n (x) = 1) = p n . 

oo / oo \ 

If E] p n < oo, then p ( E^ X n = oo ) =0. 

n=l \n=1 / 

We use the setting where X is P 1 Ai(Si i i), p is and T is the set of 

measurable sets. Let {Y n } be random variables defined by — log £ sy s(E t ( n ){ ■ )). 
The value p(Y n > — loge) stands for the measure of the set of points ([X],A) in 
T 1 M(Si,i) such that the systole length is shorter than e after earthquaking by 
time t(n) in the direction A. Using the fact that the measure v is invariant under 
earthquake flows and Proposition 01 

p(Y n > — loge) = O ( —^— ) as e —> 0. 

Vloge/ 

Let {r ra } be a sequence of positive real numbers and 




n 


1 r n , 

0 otherwise. 


It suffices to show that for all S > 0 and almost every ([X], A), Y n > (1 + <5) logn for 
only a finite number of n. This follows from setting r n = (1 + <5) logn and applying 
Theorem 0 A quick computation shows that 

E Pn = E M > (1 + S) log n) = Y, ° ( (1 + S) w n ) <0 °■ 

Hence p (X)^Li X n = oo) = 0, which is equivalent to p X n < oo) = 1, and 

this implies that the limsup of Y n / logn is bounded above by 1 almost everywhere. 

□ 


It is important to note that f(n) is a sequence chosen prior to the almost every¬ 
where property. For example, a constant sequence would satisfy the upper bound 
trivially for every ([X], A). The choice of t(n) = n were used in classical logarith¬ 
mic law papers Emn] and we will show that the same choice does not hold for the 
earthquake flow. 


5. Main Theorem 

We now prove the main theorem. The general notation is as in the proof of 
Proposition 0] 
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Theorem |TJ Let S be the once-punctured, torus or the four-punctured sphere. For 
almost every ([X],A) € V 1 M.(S) with respect to v{S), we have 


lim sup 

t—¥ OO 


-log £ sys (E t (X, A)) 

Vt 


= OO. 


Proof. We consider A with irrational slope. Recall that for each fixed e > 0, T > 0, 
and to > N(T, e), there exists a finite sequences {Xj}j=o,...,k that corresponds with 

the sequence of rational twists {Pj/<lj}j=o . k and the sequence of time {tj}j=o,...,k- 

Let a be a closed geodesic that realizes £ sys (X, A) and 7 = for some m> N 
large enough. 

Proposition[2]shows that — log £ sys (E tj ( X , A)) is approximately £ 7 (X)/(2qj). To 
finish the proof we estimate tj for almost every (X, A). For the earthquake flow 
E t {X, (£ 7 (X)) _ 1 7 ), we have 

tj ~ to = (— ~ —) (£y(X)) 2 , 

\Qj QoJ 

where t 0 , corresponding to X 0 , could be negative. The factor of (£ 7 (X)) 2 is the 
result of our normalization of the twist parameter and the weight on 7 . 

The value of to is bounded by half the maximum time that a realizes the systole 
length along the earthquake flow. A simple upper bound is obtained by considering 
the time it takes to go from po /qo to pi/qi- Again using properties of continued 
fractions from 0] and notation from the proof of Proposition^ we have the following 
bound that compares convergents. 

P± _ Po < 1 

Qj Qo Q j h(j)-\-iQj 

Combining our approximations above, we obtain 

— log £sys{Et j (X, X)) ^ _ l 7 pQ/(2<a) - £ _ 

^ ^ 7 (X))V(a M i) + iqI) + (£,(X)) 2 /(a hU)+1 q 2 ) 


We use the fact that q\ is much larger than qj to simplify into 
— \og£ sys (E t . (X, A)) 1 _ 

yj- > 2 V a ' h (j)+ 1 £ 

for j large enough. 

As to goes to infinity, the slope of 7 ™ with respect to the chosen basis converges 
to some irrational number s, the slope of A. Recall that po/qo depends on A' and 
a choice of a systole curve a. As a function of to, po/qo converges to some value 
a = cr(s). In fact, for any fixed j, the sequence {pj/qj} converges to a as to —> 00 . 
This comes from the observation that any finite sequence {p 3 {rn)/qj (to) }j = i is a 
finite subsequence of the convergents of a for m large enough. 

By choosing T arbitrarily large and e arbitrarily small, we get 


lim sup 

£—KX5 


-log £ sys {E t (X, A)) 

Vt 


> lim sup i Ja h (j) +1 (a) - e, 

j—yoo £ ’ 


where a n (a) is the partial quotients of a. The right-hand-side is infinity for almost 
every real number cr. 

The final step is to analyze the function cr : K. —► R that depends only on X and a 
choice of a systole curve a. Rational slopes would be mapped to the rational twists 
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in the sense of po/qo with respect to 7 m . To be explicit, suppose that 7™, = 7 (a, b), 
a = 7(c, d), and a(a/b ) = po/qo- The value of qo is equal to the intersection number 
of 7(0, b) and 7(c, d), that is, qo = \ad— bc\. The value of po depends on X since it 
is related to the choice of twist parameter. The function a is a continuous bijection 
over the rational numbers and irrational numbers are mapped by continuity. 

For any integer q we consider the set of integer pairs (a, b) such that \ad—bc\ = q 
for (c, d) fixed. This set is mapped to the set p/q under a. Hence the derivative of a 
is controlled by the choice of (c, d), which is based on a choice of basis for 717 (Sip). 

In an open neighborhood of a fixed (X, A), we can pick a suitable choice of basis 
for the fundamental group of S\ t i so that cr is an absolutely continuous function 
over the set of slopes for A in the neighborhood. This implies that sets of full 
measure are mapped to sets of full measure under a. Therefore for almost every 
slope s, ct(s) satisfies 

hm sup a h ( j)+ 1 (cr{s)) = 00. 

00 

Therefore by applying Fubini’s Theorem to the measure of the set of ([X],A) 
such that C^i = 00, we obtain that 


lim sup 

t—¥ OO 


-log £ sys {E t (X, A)) 

Vi 


= OO 


for almost every ([X],A) £ V 1 M{S) with respect to v(S). 


□ 


The special cases when a is equal to the golden ratio or when the tail of the 
continued fraction expansion is all l’s give us the corollary below. 


Corollary 8. For all ([X], A) e V l M{S), 

Urn sup ~<°gC„.(E.(UA)) .. 

t—> 00 y/t 2 yj (f>^ — 1 


= 


1 + V5 


if and only if A has irrational slope. Furthermore the inequality is sharp. 

Proof. If A has rational slope, then 

lim sup [ - log £ sys (E t (X, A))] 

t—¥ OO 

is a constant that depends only on ^(X). 

For irrational A we consider er(A) and its continued fraction expansion a m /b m . 
The value of C^ is larger for faster growing 6 m ’s, i.e., faster decreasing g/s. We can 
focus on the case when cr(A) is equal to the golden ratio ip since it is the irrational 
number with the slowest possible growth of b m . 

We estimate C^ similar to the steps in the proof of Theorem [ 1 ] We use two 
properties of ip. 

Pf _Po = 1 1 + 

Qj 9o qj ■ 9i-2 9j -2 • 9j-4 

where the last term is either l/go92 or l/go9i- The ratio between values of qj is 
controlled. 

bn 


lim 

m—too 0 


m+c 


= F ■ 
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Combining the approximations all together we have 


-\ogt sys {E tj {X,X)) 

\fh 


l 7 pQ/(2<fr) 



As qj’s and j become larger, we get our desired result for er(A) = ip. 


lim sup 

t—¥ OO 


-log e sys {E t (X, A)) 

Vt 


2^ 4 -l' 


□ 
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